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ABSTRACT 

We investigate how the current and future uncertainty on the Hubble constant affects 
the uncertainty in the Equation of State of Dark Energy (w) and the total density 
of the Universe (fltot)- We start with the approximate linea r relations between the 
cosmological parameters as presented by (jSpergel et a/.l 120 06). and use the standard 
error-propagation relations to estimate the effects of improving the Cosmic Microwave 
Background (CMB) parameters as well as the Hubble constant (-Ho) on °ur knowledge 
of the Equation of State (EOS) of Dark Energy. Because we do not assume a flat 
universe, we also estimate the attainable accuracy on f2 tot and the spatial curvature 
of the Universe. 

In one limiting case, we assume that the constraints provided by additional data 
(galaxy clustering, weak lensing and so forth) do not improve significantly, while the 
error on the Hubble constant is decreased by a factor up to ten. The other limiting 
case of significantly improved additional data with current Hp errors has been investi- 
gated by the Dark Energy Task Force [DETF: lAibrecht et all (|2006t )]. For the former 
scenario, we find that future improvements of the determination of the CMB hardly 
changes the accuracy with which the EOS and Qtot are know, unless the Hubble con- 
stant can be measured with an accuracy of several percent. The conclusion of the 
DETF is that the Hubble constant hardly matters if the additional data is sufficiently 
accurate. We find that a combination of moderate Hq improvements with moderately 
improved "other" data might significantly constrain the evolution of dark energy, but 
at a reduced cost. 

We review in some detail several methods that might yield extra-galactic distances 
with errors of order of several percent, where we focus on the current and future 
strengths and weaknesses of the methods. Specifically we review the following: the 
Velocity Field method, two Maser methods, four Light Echo techniques, two Binary 
Star methods, and the "Rotational Parallax" (RP) technique. Because these methods 
substantially rely on geometry rather than astrophysics or cosmology, their results are 
quite robust. 

In particular we focus on the advantages of the Rotational Parallax technique 
which can provide accurate (1%), single-step, and bias-free distances to Local Group 
galaxies. These distances can be used to improve the zero-point for other distance 
indicators which in turn would then be able to determine the Hubble constant to 
greater accuracy than they currently do. Achieving an accuracy of a few percent in 
the zero-point distances to M 31, M 33 and the LMC by the RP method requires radial 
velocities at the 10 kms -1 level and proper motions attainable by future astrometric 
missions such as SIM, GAIA and OBSS, or by future radio observatories such as the 
SKA. 

Key words: galaxies: distances and redshifts, Local Group — galaxies: individual 
(Large Magellanic Cloud, M31, M33) - cosmology: cosmological parameters, dark 
energy and distance scale - astrometry and celestial mechanics: astrometry 
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1 INTRODUCTION 

A Key Project of the Hubble Space Telescope (HST) 
resulted in a fairly accurate determination of the Hub- 
ble constant: Hp = 74 ± 2 (rando m) ± 7 (systematic) 
kms^'Mpc" 1 l|Freedman et al.ll200ll ). There are a number 
of sources of systematic error, but the largest contribution to 
the systematic error is due to the uncertainty in the adopted 
distance to the Large Magellanic Cloud (± 6.5%). Notwith- 
standing the systematic uncertainty, relative distances mea- 
surements to galaxies can now be made with a precision (not 
accuracjQ) of 5-10%. 

The drive towards accurate cosmological parameters 
is in part motivated by the connection with fundamen- 
tal physics theories that might be able to predict observ- 
ables such as the value of the critical e nergy-de n sity o f 
the Universe p fr it = 3Hi/(8nG) [e.g., ICarroll I l|200ll ): 
lAlbrecht etahl (2006) and many others]. If one assumes a 
flat ACDM model, then the fluctuations of the cosmic mi- 
crowave background (CMB) as observed by WMAP imply a 
simila r value and uncertainty of Ho [—73 ± 3. ISpergel et al\ 
<|2003l . 120061 ); h ereafter referred to as WMAP07] as the 
HST results bv iFreedman et all (|200lh . If the flatness as- 
sumption is abandoned, WMAP hardly constrains Ho be- 
cause the WMAP data constrains the actual matter den- 
sity: that is to say, the product of the normalized matter 
density (fi m ) and h 2 (WMAP07, their figure 20), where 
h = 100 kms -1 Mpc^/flo, and £l m h 2 = cu m (see Obvi- 
ously, one way to determine Q m is to accurately determine 
the Hubble constant. 

Another way to determine Q m , the Hubble constant 
and other cosmological parameters is to employ "other" 
data sets that exhibit different dependencies on h. Such 
a method is preferred while the error on h is of order a 
factor two. The constraints from other data sets used by 
WMAP07 are: large-scale structure observations galaxy red- 
shift surveys), distant type-la Supernovas, Big-Bang nucle- 
osynthesis, Sunyaev-Zel'dovich (SZ) fluctuations, Lyman- 
q forest, and gravitational lensing. However, the physical 
processes underlying these constraints can be more com- 
plex than those of the CMB and/or "geometric" determi- 
nations of the Hubble constant. Reliance on these addi- 
tional data can potentially lead to s erious biases in th e 
derived cosmolog i cal pa rameters [e.g., ISeliak et al. I (120031) : 
lEfstathiou et alA l|2005l )]. Because it is beginning to become 
possible to achieve highly accurate "geometric" determina- 
tions of Ho, it will soon be more profitable to use Ho as an 
independent constraint rather than as a to-be-fitted param- 
eter. 

Combination of WMAP and the other data sets yields 
strong evidence that the Universe is indeed close to being 
flat with a total normalized density of fltot = 0.996 ± 1.6% 



1 In this paper we make a distinction between precision and ac- 
curacy, in the sense that the former refers to small internal errors, 
while the latter indicates small external errors. 



(see section [3] below) . Combining the WMAP data with the 
Hubble constant yields the equation of state (w) of Dark 
Energy (DE): w = -0.95 ± 12% [eqn. JT3}], where the equa- 
tion of state (EOS) is the ratio of pressure (p) and density 
(p). This results follows from the assumption of a flat uni- 
verse. If flatness is not assumed, our approximation yields 
w — 0.98 ± 7.2% [eqn. 0], which is close to the rigorousl y 
derived value (w = -1.062±g;j£|) of lSpergel etali i|2006l ). 
In this paper we explore how the uncertainties on Qtot and 
w depend on the uncertainty of the Hubble constant. We 
develop our models for the case that w does not evolve, and 
that the accuracies of the other data sets do not improve. 
We then re-scale these models to take w-evolution and data- 
quality improvements into account. 

We find that accurate knowledge of these parameters 
requires better CMB data (as provided by PLANCK), but 
also significantly smaller errors on Ho (SHU). In contrast, the 
Dark Energy Task Force (DETF) finds that the errors on 
Ho hardly matter as long as the othe r data sets improve sig- 
nificantly in quality [Albrecht et~"ai1 l|2006l ). and § g53] . We 
suggest that DETF-Stage IV Dark-Energy knowledge might 
be achieved with less focus on large-scale surveys of the dis- 
tant Universe by incorporating a highly accurate measure 
of Ho obtained at earlier stages such as those provided by 
the "Water Maser Cosmology Project" iBraatz et~al\ (|2006T ) 
and/or independently by a rotational-parallax program Sj4] 
If that is indeed the case, substantial cost and time sav- 
ings might be achieved. Of course, a locally determined Ho 
value also serves as a completely independent check on "tra- 
ditional" CMB-based H determinations such as WMAP07. 
Also, future astrometric missions will provide accurate con- 
straints on the minimum age of the Universe as obtained 
from binary stars (JdJ). 

This paper is organized as follows. In section [5] we try 
to summarize the current state of the zero-point of the extra 
galactic distance scale, with an emphasis on methods that 
may be able to yield extra-galactic distances at the percent 
level (§ 32TTJ ■ ln the appendix, we work out that distance 
estimates based on water-masers may be biased if elliptical 
orbits are present in AGN disks. In section [3] we discuss our 
modeling procedure to estimate the errors on the equation 
of state of Dark Energy as a function of the errors on the 
Hubble constant. Section [3] describes the Rotational Paral- 
lax method in detail, where we also differentiate between 
the abilities between the on-going and planned astrometric 
missions: GAIA, SIM and OBSS. We discuss and summarize 
our results in section [5] Here we also describe other (future) 
astrometric constraints on cosmology. 

2 THE EXTRA-GALACTIC DISTANCE SCALE 

During the past decades significant progress has been made 
on the calibration of the extra-galactic distance scale, 
and the determination of the Hubble constant with meth- 
ods such as type la Supernovas, the Tully-Fisher rela- 
tion, surface brightness fluctuations, the Tip of the Red- 
Giant Branch (TRGB) and the Fundamental Plane. These 
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distance-determination techniques are now relatively free of 
systematic uncertainties. However, the primary calibration 
of these methods comes from the period-luminosity relation 
for Cepheid variables. 

However, the true value of zero-point of the Cepheid 
distance scale is still debated, in particular the difference 
between the zero-points based on Cepheids in the Large 
Magellanic Cloud (LMC) a n d those in the Milky W a y [e.g . , 



Feast fe Catchpolel (Il997l); IMadore fe Freedman 



Reid I (Il999l) 



Groenewegen fc Oudmaiierl ( 2000l ); iDi Benedetto I 



1998 



2002) 



ments ( ESA 



Of the Cepheids with HIPPARCOS measure- 



19971 ), all but 2 have signal-to-n oise (S/N) val- 



ues less than five l|Feast fc Catchpole Ill997h . so that sys- 
tematic biases {e.g., Lutz-Kelker correction) are important. 
Thus, the calibration of the Cepheid period-luminosity re- 
lation relies on Galactic Cepheids in open clusters, and 
is therefore tied to main-sequence fitting for clusters cal- 
ibrated t o the Hyades, and th e till-recently-problematic 
Pleiades ijSoderblom et al~ll2005l ) . However, steady progress 
is being made in the calibration of Galacti c Cepheids dis- 
tances via HST trigonometric parallaxes [Benedict et al.l 
l|2002. 20071 )] and interferometric calibration of the Baade- 
Wesselink method (lLane et al.l |2002|; iNordgren et al.l l 2002l ; 



iKervella et"aill2004al ; lKervella et al. U2004bl ) 



Because of uncertainties in the Galactic calibra- 
tion, most extra-galactic distance scale studies have been 
calibrated relative to the nearby LMC. However, the 
metallicity of the LMC is substantially below that of 
those distant galaxies (and the Milky Way [MW]) that 
are used to calibrate the Supernova la distance scale 
onto the Cepheid distance scale [e.g., ISandage et al\ 
120061 )], while there are also strong indications that the 
Period-Luminosity (-Col or ) relation is non- li near [e.g., 
iNgeow fc Kanbur l i|2006l ); iGroenewegen et al\ (|2004l ) and 
references th er ein] and dependent o n metallicit y [e.j 



Gould I (119941); ISasselov et al. I (119971); ISakai etal I 1)2004 
Groenewegen et al\ ( 20041 ); ISandage et al\ ( 20061 )]. In fact, 
the range in recently-published values of Hq (58 - 75 
kms -1 Mpc" 1 ) may be entirely attri butable to difference s 
in the applied metallicity corrections (jSandage et al.| [2006). 

Several independent methods for measuring distances 
are commonly applied to the LMC. These include Cepheids, 
the Red Clump, eclipsing binaries, SN1987A, the magnitude 
of the TRGB, RR Lyraes, and Miras. Unfortunately, the full 
range of most of the distance modulus (fiL Mc) to the LMC 
extends from 18.2 to 18.7 mag (|Reid Hl999T l. The wide range 
of moduli relative to the quoted internal errors indicates 
that systematic err ors stil l domi nate the determinations of 
the LMC distance. lAlves I (|2004l ) finds that recent applica- 
tions of a variety of methods seem to settle on \ilmc = 18.5 
± 0.1 mag (or Ho = 71 ± 10 kms - Mpc - ), which corre- 
sponds nicely with the WMAP value. However, it is possible 
that such convergence is partly related to (over-) confidence 
in the WMAP results. After all, some inconsistencies re- 
main. For example, the K-band Red-Cl ump distance to the 
Galactic centre equals 8.24 ± 0.4 kpc (I Alves I I2OO0T ) while 



geometric methods yield values one to two sigma smaller, 
or about 7.3 kpc: 1) the "expandi ng water ma ser" distance 
to Sgr B2 equals 7.2 ± 0.7 kpc (iReid I Il993l), and 2) the 
bias-free "orbital parallax method" l|Armstrong et al.lll992l ) 
of stars orbiting Sgr A* yields 7.4 ± 0.2 kpc (Ghez, 2006, 

private communications). 

The TRGB method i|Da Costa fc Armandroffl Il990l ; 
iLee et al. I Il993l ) is based on a break in the luminosity 
function of the "giant branch," which is due to to the 
maxi mum attainable luminosity before the helium flash 
[e.g jiben fc Renzini I (|l983l ): ISalaris fc Cassisi I (|l997l )]. The 
TRGB method is independent of the adopted distance to 
the LMC because it uses Galactic Globular Clusters (and 
hence RR Lyrafl) as calibrators instead. For metal-poor 
systems, this method has a metallicity dependence of ap- 
proximately Mf^ GB » -3.66 + 0.48[Fe/tf] + O.U[Fe/H} 2 
jBellazzini et al.ll200ll ) which is fairly small for [Fe/H] < - 
0.6, but increases rapidly towards higher metallicities. If the 
metallicity can be determined independent!^ the random 
errors of the TRGB method correspond to a distance error 
of several percent, while the systematic error of approxi- 
mately 10% is mostly due to the uncertainti es in the Globu- 
lar Clu st er distance scale [s ee, for example. I Mouhcine et all 
|2005bl ); iRizzi et all (|2006l ) for recent discussions]. Thus, 
the TRGB will provide accurate extra-galactic distances 
as long as a sufficien t num b er of RGB stars ] > 50 , e.g., 
IMadore fc Freedman I (|l995h ; | Mak arov et al. (2006)] can 
be resolved, and the metallicity can be estimated accu- 
rately. However, given that galaxy-halos are often the target 
for TRGB programs a nd that these halos often have large 
spread in metallicity (iDurrell et al.l l200ll ; iMouhcine et al.l 
20053; IMouhcine 1 I2006I ). metallicity effects will be impor- 
tant for very accurate TRGB results. Nevertheless, a large 
majority of accurate galaxy distances with D < 12 Mpc cur- 
rentl y derives from the TRGB method (|Karachentsev et al. I 
2006). In fact, the TRGB method may be the most robust 
photometric method available for distance determination in 
the nearby Universe. 



Although it would be preferable if the first steps 
of the distance ladder could be avoided altogether via 
trigonometric parallaxes, such measurements are beyond 
the capability of future astromet ric missions such as 
NASA's SIM PlanetQuesfl [e. g., (|Edberg et al. I 120051 )]. 
ESA's GAlffl dPerrvman ll2002l ) and the proposed OBSS 
mission l| Johnston et al.ll2006l ). The accuracy of SIMs grid 
is expected to be about 3.5 /xas, which corresponds to a dis- 
tance error of about 19% for the LMC, while for GAIA and 
OBSS the accuracy is coarser by a factor of two. Note that 



Bcllazzini ct al. I 1120011 ) employ a distance to ui Cen based on a 
double-line spectroscopic binary. 

3 In practice, the metallicity is often estimated from the photom- 
etry; 

4 http:/ /planetquest.jpl. nasa.gov/documents/WhitcPaper05ver 18_final.pdf 

5 http:/ /astro. estec.esa.nl/SA-general/Projects/GAIA/gaia. html 
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SIMs status is current l y unc leafl The Square Kilometer 
Array [SKjQ Schili zzil ri2004l)l will have astrom etric capa- 
bilities (|Fomalont fc Reid 112004 iFomalont |[2005l ) compara- 
ble to GAIA and OBSS for those sources that have strong 
radio emission, and will thus not be able to measure extra- 
galactic parallaxes directly. 

Often, many members of "standard candle" group are 
located within several kpc from the Sun. Most of these 
standard-candle methods (supernovae being an obvious ex- 
ception) will be revitalized with the advent of the next gener- 
ation space astrometry missions. Also, it is possible to obtain 
rotational parallaxes at the required acc uracy required for 
extern al ga laxies distance dete r mination [Peterson fc Shaol 
l| 19971 ) and lolling fc Peterson I l|2000l . lin preparation! ) . here- 
after referred to as OP2000]. This is discussed in §|3] below. 



2.0.1 "Nearby" Extra- Galactic Calibrators 

A future determination of an accurate distance to a single 
galaxy does not amount to an accurate determination of the 
Hubble constant. In most cases, the Hubble constant would 
need to be determined via secondary calibrators in galaxies 
at larger distances because typical deviations from the linear 
Hubble flow due to interactions with companions, groups, 
(super) clusters and voids do occur. For nearby galaxies, the 
size of these effects virtually guarantee that the measured 
recession velocity of the calibrator is atypical for the dis- 
tance of the galaxy. If a typical deviation from the Hubble 
flow is, say, 200 kms -1 , then a 1% accuracy of the Hub- 
ble constant would require a calibrator galaxy at redshift 
200/0.01 = 20,000 kms" 1 , or about 280 Mpc to eliminate 
the uncertainty induced by non-linear Hubble flow. 

Currently, the only galaxy with a distance error 
of ~5% is NGC 4258. At a distance of 7.2 Mpc 
l|Herrnstein et al. 1 1 19991 ; iHumphrevs et al. 1 120051 ) has just 
begun to serve as a calibrator for secondary distance 
indicators: the TRGB method (|Mouhcine et al.l l2005bl ) 
and the Cepheid Period-Luminosity-A bundance relation 
l|Newman et all 1200 ll ; iMacri et al.ll200rj) . However, a closer 
calibrator with at least comparable accuracy is preferable 
because secondary effects due to metallicity and so forth 
can be studied and corrected-for on a star-by-star basis. 
The LMC, M 31 and M 33 are the obvious targets, while the 
rotational parallax method is the ideal method to obtain 
unbiased distances [see § 3] . 

Alternatively, one can use a larger number of less ac- 
curate nearby calibrator galaxies. If Cepheids are used as 
secondary distance indicators, the resulting distance error 
might not decrease much below 5% due to metallicity ef- 
fects as well as uncertainties in the Helium abundance and 



6 While it is "placed on hold" there is 20 - 30 M$ (US) per year 
in NASA's budget for the SIM mission. 

7 http://www.skatelescope.org/; The SKA is planned to be fully 
operational around 2020. 



other peculiarities of individual Cepheids (I thank the ref- 
eree for pointing this out). On the other hand, it is quite 
possible that these individual peculiarities can be calibrated 
to better accuracy with highly accurate parallax and photo- 
metric GAIA data for many hundreds of Galactic Cepheids, 
TRGB stars, RR Lyrae and so forth. 

2.1 "Direct" Distance Measures 

Fortunately, other more direct methods that can yield reli- 
able distances have been proposed (and applied) in the lit- 
erature. These methods can be applied in either the optical 
and/or the radio regime. Below we review the most promis- 
ing of the candidates that we have been able to identify in 
the literature. Some science projects require extra-galactic 
distances with an accuracy of 5% or 10% that will be readily 
attainable for a large number of galaxies may be sufficient, 
while other projects require the very best possible accuracy. 
The TRGB method might be applied for projects of the 
former type as soon as its zero-point has been accurately 
established. The review below is focussed on the latter type: 
methods that can deliver distance accuracies at the percent 
level. 

2.1.1 Velocity Field Methods 

The "Rotational Paralla x" (RP) method was conceived by 
iPeterson fc Shaol Jl997) in the SIM context as an extra- 
galactic variant of the "orbital parallax" method. In the 
latter method, radial velocities and proper motions of a re- 
solved binary system are combined to yie ld all the orbital 
param e ters, as well as the distance [e.g., . lArmstrong et all 
l|l992h : iDavis et al.l (|2005l )]. In the simplest extra-galactic 
application, one assumes that: 1) the orbits are close to cir- 
cular, 2) with deviations only due to (small) random mo- 
tions, 3) the inclination can be obtained from the analysis 
of the radial velocity field, and 4) the dynamical centre is 
known from the velocity field. With these simplifying as- 
sumptions and the use of just two stars per galaxy, the sys- 
temic motion of the galaxy can be subtracted and the dis- 
tance determined from the orbital model and the observed 
radial an d proper motions. A more detailed treatme nt can be 
found in l|011ing fc Peterson Il2000l lin preparation! ) and § |U 
The trigonometric-, orbital- and rotational parallax meth- 
ods rely only on geometric relations, Newton's laws of mo- 
tion (and the ability to measure accurate positions, proper 
motions and radial velocities). These methods are entirely 
independent of photometric transformations, metallicities, 
effective-temperature scales, pulsation theory and other un- 
certain prope rties of stars such as their interiors and their 
atmospheres l|Kuruczll2002l ). 

The "Velocity Gradient" method is in part geometric, 
and in part photometric and has been appli ed to op e n clus - 
ters and modified to work for the LMC [see iGould I |2000), 
and references therein]. In this method, the transverse ve- 
locity of the centre-of-mass of the cluster/galaxy induces 
an angular difference between the kinematic and photomet- 
ric major axes (about 20° for the LMC). In practice this 
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method would require radial velocity and proper motion 
measurements of about 10,000 stars with accuracies better 
than about 10 kms -1 and 150 /iasyr -1 , respectively, which 
are easily attainable by GAIA. The most difficult part of this 
method is the determination of the photometric major-axis: 
this requires 0.1 mmag photometric accuracies over roughly 
10 de grees as well as accurate reddening estimates (I Gould I 
2000). Although it may not be possible to obtain such data 
from the ground, GAIA is well-equipped to deliver all the 
required data for this method. This method might also be 
applied to GAIA data of M 33. The high inclination of M 31 
will likely preclude successful application of the velocity gra- 
dient method. Note that in practice the Rotational Parallax 
Method (or a simplified version thereof) may be needed to 
det ermine the systemic prope r motion of the LMC [see § [4] 
and lOlling fe Peterson I (|2000ri ] . 

2.1.2 Maser Methods 

Water masers associated with star forming regions scattered 
across the disk of a spiral galaxy can be used to estimate 
the distance to that galaxy. With two masers sources lo- 
cated on "opposite" sides of the centre, the known rotation 
curve and in clination of the ga l axy c an be used to derive 
the distance. iBrunthaler et alTI l|2005l ) identified 37 maser 
features in two masing complexes in M 33 for which they 
could determine proper motions. From the per-complex av- 
erage motion in c o mbina tion with the velocity-field model, 
IBrunthaler et al. I (2005) determined the distance to M 33 
with an error of 23%. Note that this method is a version of 
the rotational parallax meth od, but one that relies on exter- 
nal information [see § HI and lOlling fc Peterson! (|2000T) ] . All 
the systematic effects that could affect the RP method are 
also important for the water-maser work, especially when 
only a limited number of sources are available per galaxy 
[M 31 i s expected to have 17 H2O masers at SKA sen- 
sitivity l|Brunthaler et al.l [20061 1]. When the SKA becomes 
operational, many more fainter m aser sources will be de- 
tected in galaxies out to 100 Mpc jFomalont II2005T I. while 
a distance accuracy of 1 % can be achieved out to 30 Mpc 
ijFomalont fc Reidl [2004'). At this distance, the perturba- 
tions in the Hubble flow would allow for a determination of 
about 10% per galaxy: the more distant sources would allow 
for the 1% goal on Hq. 

Rather than averaging the proper motions per masing 
region, the H2O maser data can be utilized in a manner sim- 
ilar to their Galactic equivalents to derive distances. That 
is to say: a dynamical model of the star-forming region is 
used to predict the relation between the observable radial 
velocities and proper motions. Such models can either be 
based on ordered motions (such as rotation or expansi on) or 
on isotropic random motions. See lArgon et alj ([2004) for a 
full description, references and an application to the IC 133 
star- forming complex in M 33. For the case of th e 14 year 
observing campaign of IC 133, lArgon et all l|2004[ ) adopt a 
random+systematic error of 22%, which is dominated by: 1) 
the unknown centre of expansion (±16%), 2) random error 



(±12%), and 3) an assumed expansion velocity (and error) of 
the group of maser spots (as determined by comparison with 
Galactic H2O masers; ±9%) and an assumed systemic ve- 
locity (and error) based on co-incident CO emission (±4%). 
Although the expended effort is already enormous, it seems 
likely that the results could improve by at least a factor 
of two as a result of continued long-term monitoring (espe- 
cially with the SKA), which would reduce the RMS errors 
and would eliminate the secondary possibility for the expan- 
sion centre. 

NGC 4258 is the current record holder for the most 
accurate extra-galactic distance, with a total error of 
about 5% as derived from the ge ometry and dynamics of 
its nuclear water mase r sources (|Herrnstein et al~l 1 19991 : 
iHumphrevs et al~l 12005). The H2O masers are thought to 
arise in an almost edge-on circum-nuclear disk in Keplerian 
rotation. In this model, the masers occur close to the lo- 
cations where the line-of-sight velocity gradients are small: 
that is to say, close to the major and minor axes of the 
velocity field (the major axis is also called the "midline"). 
Maser spots close to the major and minor axes are called the 
"high-velocity" and "systemic" masers, respectively. With 
the assumption of circular orbits, the observed positions and 
velocities can be converted to a model where the inclination 
(i), the position angle (<j>) an d the rotation speed (V c ) de- 
pend on radius (R) only [see iHerrnstein et all (|2005l ) and 
references therein]. Given this model, the distance can be 
derived in two independent ways from the systemic masers. 
First, the circular speeds of the systemic masers are directly 
observable as tangential proper motion with a value of about 
30 ^tasyr -1 . Because phase-referencing techniques are used 
to obta in the VLBA positions, absolute position information 
is lost (|Herrnstein et afl feoOS). so that the systemic proper 
motion of the galaxy is effectively subtracted from the ob- 
servations. Thus, the above relations amount to a special- 
ized version of the rotational parallax technique (see Q. 
The conversion from observed position to radial coordinate 
R introduces the largest systematic uncertainty in distance 
because V c , i and cj> change fairly rapidly with radius. The 
second technique employs the change in radial velocity (V r ) 
that occurs due to the orbital motion of the systemic masers. 
In almost edge-on geometry and with V c ~ 1100 kms -1 and 
P about 750 yr, the accelerations are about 9 kms -1 yr _1 . 
These two geometric techniques yield consistent results. 

2.1.3 Nuclear Water Maser: Eccentric Orbits 

A systematic uncertainty that has hitherto not been inves- 
tigated is that of elliptical orbits in the masing region. It 
is well established that the majority of the high-luminosity 
cousins of NGC 4258 often exh ibit significant degr ees of 
non-axisymmetry in their disks (IStrateva et al.| [2003). Fur- 
thermore, models employing disks with substantial elliptic- 
ity were very successful in describi ng the observed optica l 
double-peaked line profiles of AGN l|Eracleous et al.iri995l "). 
who infer orbital eccentricities of order three-tenths. Such 
large eccentricities could result if the AGN contains in fact 
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a binary black hole, or the tidal disruption of a normal star 
for the case of a single black hole. In such models, it seems 
likely that the water-masing region would also be affected 
because the radio emission is generated just a few times 
further out than the optical emission lines. However, other 
explanations exist for the optical emission lines, and it is 
not cl ear that elliptical orbits do in fact occur in accretion 
disks (|Storchi-Bergmann et a.1.1 12003T) ■ In fact, models with 
such large and persistent eccentricities have been ruled out 
on the basis of long-term spectroscopic mon itoring of several 
AGN (|Gezari. Halpern fc Eracleous 1120071 ). Nevertheless, it 
seems that the possibility of elliptical orbits should be taken 
into account in the systematic-error budget of nuclear wa- 
ter masers, because e ccentricities thirty times small e r than 
those ruled out by iGezari. Halpern fc Eracleous I (|2007l ) 
would still bias the inferred distances (see below). Observa- 
tions of other maser lines (if present) which probe different 
regions of the accretion disk would also minimize the uncer- 
tainty associated with large-scale non-circular motions. 

Here we analyze a simple version of the general case of 
elliptical orbits where the major axis of the ellipse is aligned 
with the line of sight and show that this geometry yield 
observables that are indistinguishable from the case of cir- 
cular orbits, albeit for a different distance, inclination and 
mass. For illustrative purposes we also assume that the disk 
is not warped. This case is worked out in some detail in 
Appendix [X] where we find the following relation the dis- 
tance [eqn. [fA20])]: D c /D t = ^/(l q= e) 3 / (1 ± e) « 1 =p 2e, 
where D c is the distance that would be inferred under the as- 
sumption of circular orbits, Dt the true distance as derived 
for non-circular orbits, and where the minus (plus) signs 
correspond to case that the masing occurs at peri- (apo-) 
centre. Thus, if circular orbits are assumed but the orbits 
have in fact a small ellipticity, the true distance is smaller 
(larger) than the inferred distance by approximately twice 
the true orbital eccentricity if the systemic masers occur at 
peri- (apo-) centre. We have not investigated the case that 
the major axis is not aligned with the line of sight, but we 
surmise that a case will result intermediate between the peri- 
and apo-centre situations, so that the true distance will be 
uncertain by a factor of approximately (1 ± 2e). 

2.1.4 Light Echo Methods 

"Light Echo" methods have been proposed and applied 
to the absorption of the super nova (SN) "pulse" by gas 
surrounding the progenitor, which is then re-emitted in 
the form of line radiation. This method has been ap- 
plied to SN 1 987A in the LMC by several author s [e.f 
Panagia et ail dl99ll ); [Gould I (|l995l ): lGould fc Uzal (|l998 
Panagia 1 19991 )]. but it leads to a sys tematic (unre solvable) 
distance uncertainty of about 10% (|Gould 1 12000| ). So far, 
this method has only been applied to SN 1987A. 

Alternatively, one can attempt to use the polarized re- 
flection of the SN pulse off dust in the curcum-stellar ma- 
terial. The polarization condition strongly favors scatter- 
ing by 90°, so that the physical radius equals the speed of 



light times the time-lag between the SN pulse and the de- 
tection of the polarized scattering event (|Sparkslll994l ). If 
the proper motion of the light echo can be determined, the 
depend ence on the a ssumed scattering angle is largely elim- 
inated (Sparks 1996). Although previous search es for such 
events have not produced th e hoped-for results (|Boffi et al.l 
ll999l ; lRomaniello et aLl [2005), current workers indicate that 
the method works quite well (Sugerman, 2006, private com- 
munications) . 

A vari ant of the "pulse " meth od has recently been de- 
veloped by iDraine fc Bond! l|2004l) . Their method employs 
strong X-ray sources behind nearby galaxies. The dust in 
the foreground galaxy scatters the X-rays to produce a halo 
around the source position. If the background source is vari- 
able, then intensity of this halo will vary with time and 
position, where the time-lag (At) is directly related to the 
angular distance (9) from the source and distance (D) to the 
foreground galaxy: At ~ 140 D 9\ 00 d ays, with D in Mpc and 
0ioo in units of 100 arcsec. For M 3l lDraine fc Bond! (|2004T l 
estimate that a Chandra observing campaign of 5C 3.76 with 
a duration of 2 months (4 months) would yield a distance 
accuracy of 4% (1%). However, less observing time would 
be required if observing starts after an "outburst" of the 
background source. This method might also be applied if a 
gamma-ray burst would occur behind a nearby galaxy. 

The "Expanding Photosphere" method can also be used 
to model the light curve and expansion history of super- 
novas. In this method the distance follows from the mea- 
sured angular size and the phy sical size calculated fr om in- 
tegration of the radial velocities. iMitchell et all (|2002l ) devel- 
oped the "Spectral-fitting Expanding Atmosphere Method" 
which is specifically tailored to the time-varying and non- 
LTE nature of SN atmospheres (ejecta). This method has 
distance accuracies of order 10% and has been applied out 
to z=0.3 (|Nugent et al.ll2006l ). 



2.1.5 Binary Star Methods 

Detached Eclipsing Binaries (DEBs) are expected to be good 
distance indi cators throughout the Local Group, possibly a s 



good as 1% (iPaczvriski fc Sasselov"1ll997l ; iPaczvhski Ill997l ; 
IFitzpatrick et al.ll2003l ; IWilson II2004I ). In this method, pho- 
tometric observations of the eclipses determine the period, 
the flux ratio and the ratio of the stellar radii (R*) and 
the size of the orbit, while radial velocity observations de- 
termine the latter quantity directly. Spectroscopy and/or 
multi-band photometry can yield the effective temperature 
(T e ff) of the stars. The absolute luminosity (L*) can be de- 
termined if the relation between surface brightness (F*) and 
T e ff is reliably determined, so that: L» = iirR 2 F*. After 
correction for extinction, the observed brightness (F i, s ) then 
yields the distance: d — y/L*/(4irF b s ) = R, y/Ft,/F bs- A f- 
ter 10 years o f active r esearch [e.g.. iGuinan et al. I d 19981): 
iKaluzny et all (Il998l): IStanek et all (Il998l); Ribas et al I 
d2002l): IFitzpatrick et alj il2003l); iRibas et all l|2004 120051 ); 
iHilditch et all (|2005l ): lBonanos et all l|2006l )]7the best accu- 
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racies are currently about 5%. Possible reasons for the accu- 
racy gap are the difficulty to determine the effective temper- 
ature, surface brightness and extinction. Also, the early-type 
main-sequence stars that are used in the DEB method have 
no local analogs with accurate distances and surface bright- 
ness measurements: none of the 150 HIPPA RCOS stars ear- 
lier than B3 have a parallax better than 5%. We expect 
that the GAIA data will provide the necessary calibration 
data (distances and extinction) to fully exploit the distance- 
scale potential of DEBs. Finally, DEBs are not ideal dis- 
tance indicators: when using semi-detached systems or over- 
contact systems there are fewer (stellar) parameters to be 
solved for so that more accurate distanc es can be obtained 
[Wilson, 2002, private communications; IWvithe fc Wilson"! 
l|2002l ): IWilsonTlj2004l )]. 

A fundamental-physics method involving binary stars is 
via the detection of gravitational waves of close white dwarf 
binary (C WDB) stars with LIS A. After 10 years of LISA 
operation. ICoorav fc Setol l|2005l ) expect to be able to obtain 
a distance error for the LMC of order 4.5% x ^JN C wdb/22, 
where the expected number of CWDBs (NcwDB-,ex P = 22) 
may be uncertain by a factor of ten. 



3 IMPLICATIONS FOR COSMOLOGY 

There are at least four reasons to determine Ho with an 
accuracy of several percent or better: 1) to determine dis- 
tances to external galaxies via redshift measurements, 2) to 
calibrate accurately other distance indicators, 3) to deter- 
mine the equation of state of Dark Energy, and 4) to deter- 
mine the total density (fltot) of the universe. The equation of 
state (EOS) of Dark Energy (DE) tells us something about 
its nature. 

Several lines of evidence suggest that we live in a Uni- 
verse that is close to critical density (fitot ~ 1). For exam- 
ple, the WMAP data and the HST constraint on the Hub- 
ble constant (h = 0.74 ± 0.08) results in Q tot ~ 0.996 and 
a mass density (baryons & dark matter) of 23%, so that 
about 76% of the energy- density of the universe is unrelated 
to gravitating matter. The most obvious candidate for this 
Dark Energy is a cosmological constant A, which is allowed 
for by General Relativity. Other forms of DE are suggested 
by theoretical physics. One difference between these various 
possibilities is the equation of state of the proposed DE can- 
didates. For example, the cosmological constant would have 
w = —1, cosmic strings have w — —1/3, domain walls have 
w = —2/3, while Quintessence can come in multiple vari- 
eties: fixed with w > — 0.8, or with a time-variable w value 
jPeebles fc Ratra Il2003h . We discuss two scenarios. In iji l3.ll 
and i) ii3.2l we discuss the case of a constant w, while in § i]3.5l 
we discuss the case of varying w, and compare our results 
with the recent concl usions of the Dark Energy Task Force 
l|Albrecht et alj|200fj) . 

Hu (2005; hereafter H05) analyzes in quite some de- 
tail how measurements of Ho and its variation with redshift 



would be beneficial for our understanding of dark energy and 
its evolution. He concludes that: "... the Hubble constant is 
the single most useful complement to CMB parameters for 
dark energy studies ... [if Ho is] ... accurate to the percent 
level ... ." However, Hu does not explicitly describes how 
errors on the dark-energy parameters are related to errors 
on Ho, while such information is of pre-eminent importance 
when justifying and planning observational campaigns to 
determine Ho, especially if space missions are considered. 
In the remainder of this section, we present some analyti- 
cal estimates for the relationship between the errors on Ho 
and the dark-energy parameters. These relations are approx- 
imate and are intended to facilitate estimating the required 
observing time (and cost) of observing programs aimed at 
determining Ho- 

For the reader's convenience, we here present a short 
summary of Hu's (2005) description as to how CMB data 
physically relate to the Hubble constant. In fact, the CMB 
experiments determine the photon-to-baryon ratio (i?*) and 
the matter-to-radiation ratio (r*) from the locations and 
ratios of the acoustic peaks. H05 finds: r* oc °^ 126 , 

and R, oc 0223 l+'z > where w 6 = Qbh 2 is the physical 
baryon density. These relations allow us to derive the er- 
rors (e) on R* and r* from the errors on u m and ut as given 
by WMAP07. The results are: e r „=7.1% and e R ,=3.6% for 
uj m = 0.126 ± 0.009, and 100 u b = 2.23 ± 0.08. Physically, 
the co-moving size of acoustic oscillations (s*) equals the 
distance a sound wave has traveled during the time between 
the Big Bang and recombination, and depends on R t , r, , 
2* and the sound speed (and hence the matter density uj m )'- 
s, « 139.8 (offi) -0 ' 252 (oms) ' 083 M P C (scaled from H05), 
with a 1% uncertainty. Note that z t depends only weakly 
on cosmology. The CMB data also yields the location of the 
acoustic peak (£a) with smaller error than that of s„. Thus, 
the angular-diameter distance of the acoustic peak (Da) is 
well-established observationally: Da = £as,/tt, while theo- 
retically Da depends on cosmology: 

D A (n K =0) = o.£^da (1) 

where we neglect the small contribution from the density in 
relativistic particles (fi^/a 4 ), and where a — 1/(1 + z) is the 
scale factor of the universe. If the curvature densi ty (Qk) 
does n ot vanish, eqn. ([1]) needs to be modified [e.g.. ICarroll I 
|200ll )]. In this simple model, the sum of the f2-terms equals 
unity so that: 

fi A = 1 - n m = 1 - ujm/h 2 w 0.770 ± 0.048 , (3) 

where the approximate value for eqn. <(3j follows from the 
CMB value for uj m and the HST value for h. Thus, accu- 
rate determinations of uj m and h would significantly con- 
strain the dark-energy. However, because the measured Da 
is an integral constraint, it reveals nothing about the vari- 
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ation of H with scale-factor (or redshift or time). If in re- 
ality w varies with scale factor, while it is assumed to be 
constant, then application of the integral constraint would 
result in an erroneous value of Ho and/or w. However, 
if w is constant, as would be the case if the DE results 
from the cosmological constant, cosmic string, domain walls, 
etc., then a determination of Ho would nail down w. Be- 
cause the variation of w is an important discriminator be- 
tween various po s sible dark-energy candid a tes [see, e.g., 
iPeebles fc Ratra I rt20O3h; iFlambauml (|2006h : lAvelino et ail 
(|2006t) ; lAlbrecht et al l (|2006l 'l. and references therein] one 
would want to determine the variation of w accurately. 

3.1 A Flat Universe 

Here we consider the case that the universe is flat, and if 
the equation of state is constant. For this simple case, we 
can determine the EOS of dark energy from the CBM data 
and the Hubble constant only. We use this example to illus- 
trate that dark-energy information is contained in accurate 
knowledge of Ho. For the case of a flat universe (Qk = 0) 
and with Sl„ = 0, eqn. ((2)1 can be simplified to read: 



H(a) 
100 



h 2 



,3(1+™) 



(4) 



If the dark energy is in the form of a cosmological constant 
with w = — 1, then the only unknown left in eqn. Q is the 
Hubble constant which can be found from the constraint 
provided by the observed value of Da [eqn. (JT}]. From h, 
fl m and J1a follow trivially. This is the basis for the reported 
determination of the Hubble constant from the WMAP data,. 
However, the inferred value of the Hubble constant (and 
hence fi m and J1a) depends on the choice of the EOS of 
the dark energy [via the l/a 3 ' 1+m '-term in eqn. @]. This 
relationship is illustrated in Fig. 15 of WMAP07, from which 
we derive: 



1.59 - 2.56 0, 



1.59-2.56 



h 2 



(•») 



where the second equality is re-written in terms of the ob- 
servables. The errors on the EOS and the density of dark 
energy for the case of a flat universe are thus: 



3% 



2 W(S;) + (it) 



(6) 



(7) 



With uncertainties of 7.1% for w m and 9.8% for h (includ- 
ing the systematic error term) the uncertainties on w and 
Qa are almost evenly divided between the two error terms. 
However, in the future this will no longer be the case be- 
cause the PLANCK missio n is expected to reduce the error 
on uj m by a factor of eight (jEfstathiou et 0/71120051 ). This in- 
dicates that in the post-PLANCK era, our knowledge of the 
value of dark energy will be fully dominated by the uncer- 
tainty on Ho if the (systematic) error on Ho is not decreased 



signifi cantly, and if no other data is considered. iBlake et al. I 
(2004) reach a similar conclusion albeit that they argue that 
the CMB data imparts a correlation between and ui m that 
will yield the same s* so that it is impossible to go below 
e w ~ 0.1 employing Ho and CMB data alone. We note that 
this assertion is at odds with H05's who concludes that u>b 
and u> m are virtually independently determined from the 
CMD data. 



3.2 Including Curvature 

When adding other relevant data sets, the various cosmo- 
logical parameters are determine better than with WMAP 
data alone. Inspecting Figure 21, Table 11 and Figure 17 
of WMAP07, we derive the following x-by-eye relations for 
f^A, O.K and w, respectively: 



= fflKA + &if a£2a 
W = CLmK + b w K^K 



(8) 
(9) 
(10) 



We estimate the coefficients to equal: BA m ~ 0.944 ± 0.011, 
&Am ~ -0.775 and oka ~ -0.0992± 0.009, b KA ~ 0.1199± 
0.0124 (the .KA-terms are determined excluding the Ho con- 
straint in Table 11 of WMAP07) and a wK ~ -0.910±0.063, 
b W K ~ 6. In order to investigate the effects of the Hubble 
constant, we combine equations (JSJ through (|10|) in terms 
of h and arrive at: 



Ov = (0.944 ± 0.011) - 0.775- 
w 0.766 ± 0.020 (4.9%) 



h 2 



O-wK + dKAbwK + ahm.bKh.bwK + bhmbKhbwK 



h 2 



(-0.826 ± 0.109) - (0.557 ± 0.058) 
-0.95 ± 0.11 (11.6%), 



(11) 

(12) 
(13) 



where the relations that are not numbered employ the values 
and errors on ui m from the CMB and h from HST. Compar- 
ing the relation for w for a flat universe [eqn. (0] with the 
general case [eqn. (|13|l ] we notice two things. First, the error 
is larger for the general case, presumably due to the extra 
degree of freedom that arises from fitting for Qk- Second, 
the dependency on the value of Ho has decreased by a fac- 
tor of almost five as evidenced by the smaller coefficient of 
the term containing h. This is probably due to the fact that 
the relations for the general case were derived employing 
the other data sets. Also note that Qj^ at equals Qff enerai to 
within the errors, and that the derived errors are approxi- 
mately the same for the two cases. 

For the general case, we can now also obtain an expres- 
sion for the total density: 



QnoK — 



f^A + O m = BAm + (&Am + 1) -pT 

h 2 



= (0.9438 ± 0.0114) +0.225 
w 0.996 ± 0.016 (1.6%) . 



(14) 

(15) 
(16) 
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Thus, without the assumption of a fiat Universe, current 
data allows for the determination of the total density of 
the universe to plus or minus 1.6%, while the EOS of Dark 
Energy is known to about 12%. The actual values of fltot 
and w strongly hint at a flat Universe with the cosmological 
constant being the Dark Energy. 



3.3 Ho and Dark Energy 

Future CMB data will be more accurate than at present: the 
eight-year WMAP data will reduce the current errors by a 
factor of two (Spergel, 2006, private communications), while 
PLANCK data is ex pected to be eight times better than the 
3-year WMAP data (|Efstathiou et al. Il2005h . 

We can estimate the effects of more accurate CMD data, 
as well as an accurate value of the Hubble constant on our 
knowledge of the EOS of Dark Energy by employing equa- 
tion (|12[) above. In order to do so, we assume that the errors 
on the a,i and bi parameters improve while the values them- 
selves do not change with improved CMB accuracy. The 
advantage of this method is that it is easily implemented. A 
full, WMAP07-style, calculation could be performed, but is 
beyond the scope of this paper. 

As discussed in fj3] the physics of the CMB implies 
that only a relation can be found between Qa and fl m - 
This leads to a very elongated confidence region. Other data 
sets are required to break this degeneracy. For example, 
the magnitude-redshift relation for Supernova la is popu- 
lar because its confidence region is almost perpendicular 
to that of the CMB data [for r eview s, see among others , 
Carroll I (l200ll);lLewis fc Bridle I j20oll;lRiess et al. 1 1|2004) ; 



Efstathiou et al. I (|2005l V IPerlmutter I (|2005l ')] . In that case, 



the SN-Ia data "selects" part of the confidence region gen- 
erated by the CMB data to determine values for Qa and 
Sl m . However, it will hardly decrease the error in Oa (for a 
given fl m ) because the SN-Ia confidence region is very elon- 
gated in the Qa direction. Thus, orthogonal constraints are 
the way to determine parameter values, while the shape and 
direction of the joint confidence region is set by the data 
set with the smallest errors. In the remainder of this paper 
we assume that the CMB measurements will be the most 
accurate data set, so that the relations between the cosmo- 
logical parameters given by equations © through (|12p and 
(|14p are approximately valid. We will also investigate the 
consequences of relaxing this assumption. 

The results are presented in figure [1] where the top 
panel shows e m as a function of improvement of our knowl- 
edge of the CMB parameters, with respect to the WMAP 
3-year data. The four lines are computed for Hubble con- 
stants with errors of 10.8%, 5.4%, 2.7% and 1.1%, from top 
to bottom. These errors correspond to improvements with 
respect to the current error by factors of 1, 2, 4 and 10, re- 
spectively. The corresponding errors on w that are attainable 
with PLANCK-like data are, 8.9%, 4.8%, 3.0% and 2.3%, re- 
spectively. These errors are computed via the standard error 
propagation relations for equation (|12|l . The part of the error 



on w that can be attributed to uncertainty in Ho (e w ,h) de- 
pends on the relative accuracy of the CMB parameters and 
Ho. Currently, the error on Ho hardly affects our knowledge 
of w. However, with decreasing errors on the CMB param- 
eters, our ignorance of the Hubble constant becomes the 
dominant contribution to the total error on w (tw,tot). This 
is illustrated in the bottom panel of figure [T] which shows 
e w ,h/e w ,tot as a function of the accuracy of CMB data and 
the Hubble constant. 

From figure [1] we infer that at PLANCK accuracy, the 
errors on w will have only slightly decreased with respect 
to the current value of 11.5%. However, the accuracy on 
the EOS of Dark Energy improves significantly when the 
error on Ho is decreased significantly. In section 13.51 below 
we estimate the effects of significant improvements in the 
quality of the other data sets and show that in those cases 
the effects of improved knowledge of Ho diminishes. 



3.4 Ho and Q to t 

The deviation from spatial flatness [(1 — Qtot) = Qk] follows 
from Equation dH]): Q K ~ 0.004 ± 0.016. Thus, the spatial 
curvature is currently known to within a factor of four. We 
applied the same error-propagation methods for Qtot as for 
w and find a behavior for the error on Qtot very similar 
to e w , but at approximately ten times lower level. Thus, a 
determination of the Hubble constant at the 1% level would 
decrease the error on the spatial curvature by a factor 16: 
from a factor four to 25%. 



3.5 Comparison with the results from the Dark 
Energy Task Force 

In their recent r eport , the Dark Energy Task Force [DETF; 
lAlbrecht et al.l (|2006l )] presented recommendations as to 
how to increase our knowledge of dark energy, and in partic- 
ular its evolution. As a starting point, they used the CMD 
data as expected to be delivered by PLANCK and added to 
that the expected improvements of a large variety of other 
proposed space-based and/or ground based investigations 
based on Baryon Acoustic Oscillations, Galaxy Clusters, Su- 
pernovae and Weak Lensing. The uncertainty on Ho used by 
the DETF corresponds to the current systematic Hp error 
as derived from the HST data (|Freedman et alJboOlT) . 

Our work presented above is based on the case of 
a fixed EOS of the dark energy, where we also assumed 
that any future increase in accuracy would be dominated 
by the PLANCK results, and that the accuracies of the 
other data sets would remain the same. In DETF-speak, 
we use "Stage I" accuracies for the other data sets and find 
e m =8.9%. From the models tabulated by the DETF, we in- 
fer e w values of approximately 3.6%, 2.4% and 1.45% for 
the combined Stage II, Stage III and Stage IV data, respec- 
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Figure 1. Top Panel: The accuracy with which the EOS 
of Dark Energy can be determined as a function of the ac- 
curacy of the CMB data (abscissa), and the accuracy in Ho 
(curved lines). The vertical lines correspond to accuracies of the 
WMAP-l-year data (left), WMAPSyr data, WMAPSyr data 
and PLANCK (right). With PLANCK-like CMB parameters and 
_ffo accuracies of, from top to bottom, 10.8% (8 kms - 1 Mpc -1 ), 
5.4% (4 km s -1 Mpc -1 ), 2.7% (2 kms-'Mpc" 1 ) and 1.1% (0.8 
kms -1 Mpc -1 ), the achievable errors on w are: 8.9%, 4.8%, 3.0% 
and 2.3%, respectively. The "error bars" represent the case that 
the slopes bi in eqns. JS} - l|10|l are randomly varied by ±25%. 
The lines in the bottom panel show the contribution of the uncer- 
tainty on Ho to the total error on w for ej^ =8/l, 8/2, 8/4 and 
8/10 kms - 1 Mpc -1 , from top to bottom. 

tiveljQ. But note that these accuracies are only reached if all 



8 The DETF considered explicitly DE models with varying val- 
ues of the EOS: w(a) = too + (1 — a)w a , with wo the value at 
the current epoch and w a the slope of the EOS versus scale- 
factor relation. For any fitted linear relation, there is point (the 
"pivot point") where the total error is minimal. In this case, w p = 
wo + (i- — a p ) w a , and a p scale factor of the pivot point. The DETF 
used a figure of merit (FoM) to evaluate the effectiveness of the 
other data sets. The FoM is essentially FoM = 1/ (tw p e Wa ), ■ For 
models such as ours that do not include EOS evolution (w a = 0), 
the "pivot point" lies at a p = 1. We thus compare our e w with 
e„ as tabulated by the DETF. 



programs per stage are combined, and that all the programs 
in the previous stages are also executed. 

Unfortunately, we can not directly compare our results 
with those of the DETF, but we can crudely scale our 
methodology to their results in the following manner. In our 
analysis, the strongest correlation (degeneracy) between pa- 
rameters occurs for w and Q.k , with a coefficient of b W K ~ 6 
[cf. eqn. ((TO])] . The effect of adding more accurate data is 
to break the degeneracies between the cosmological param- 
eters, which we implement by decreasing the value of b m K- 
After some experimentation, we find that we recover the 
DETF e w values if we decrease b m K by factors of 2.51, 3.90 
and 7.51 for Stages II, III and IV, respectively. The DETF 
tested the effects of smaller errors on Ho for Stage IV data, 
and conclude that a factor of two improvement would lead to 
a reduction of e w by at most 50%, while we find an improve- 
ment of 26%, which is consistent with the DETF results. 

Employing this calibration, we estimate the effects of 
smaller eh values as a function of the accuracy of the other 
data. The results for Stages I, II, III and IV are presented 
in the top through bottom panels of Figure [2] respectively. 
Note that the top panel corresponds to our results discussed 
above (and Fig. [1} , while the bottom panel corresponds to 
the best-case studied by the DETF. As expected, the effects 
of an improved error for Ho has less and less of an effect on 
e w as the other data becomes more and more accurate. 

Figure [2] suggests that some improvement in the other 
data sets can be traded for an improved error on Ho- That 
is to say, Stage I with e^/20, Stage II with e/i/4, Stage III 
with en/2 and Stage IV with eh/1 yield similar accuracies 
for the equation of state of dark energy. [We give the exact 
numerical results in the caption of Figure [5]] The DETF 
quotes that each of the Stage III projects would cost tens of 
millions of dollars (US), while the Stage IV projects might 
be as expensive as 300 - 1,000 M$ each, or a total of at 
least 2,000 M$. However, the primary motivations for each 
of these Stage IV projects are not necessarily related to 
dark-energy science. Nevertheless, we suggest that signifi- 
cant cost and time savings might be obtained if more effort 
is expended on determining a more accurate Ho, or at least 
that it would be worthwhile to investigate the matter more 
thoroughly. For example, a "VLBA" effort to find and an- 
alyze more extra-galactic water masers may reduce eh by a 
factor of two to five at a cost of maybe several M$, while a 
SIM-based rotational parallax program can reduce eu by a 
factor of 5-10 and might co st 20 M$. In fact, su ch a "VLBA" 
effort is already underway (jBraatz el al . 2006). 



4 ROTATIONAL PARALLAX DISTANCES 

The method of "Rotational Parallaxes" (RP) combines 
proper motions and radial velocities of stars in external 
galaxies to yield bias-free single-step distances, analogous 
to the orbital parallax technique. The attainable accura- 
cies, of the order of one percent, are well matched to the 
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Figure 2. In each panel we present the accuracy with which 
the EOS of Dark Energy can be determined as a function of 
the accuracy of the CMB data (abscissa), and the accuracy in 
Ho (curved lines). The vertical lines correspond to accuracies of 
the WMAP- 1-year data (left), WMAPSyr data, WMAPSyr data 
and PLANCK (right). The difference between the four panels 
is the accuracy with which the other data sets are determined. 
From top to bottom, these accuracies correspond approximately 
to stages I, II, III and IV as defined by the Dark Energy Task 
Force (see text, ^3. 51 1. respectively. In each of the panels, the four 
lines correspond to the same errors on Ho: ejj =8/l, 8/2, 8/4 and 
8/10 kms" 1 Mpc -1 , from top to bottom. The attainable errors 
with PLANCK data are: e w = 8.90%, 4.81%, 3.02% and 2.27% 
for Stage I (as in Fig. [TJ e w = 3.60%, 2.06%, 1.43% and 1.20% 
for Stage II; e w = 2.41%, 1.48%, 1.14% and 1.02% for Stage III; 
6,0=1.45%, 1.07%, 0.95% and 0.91% for Stage IV. 



requirement that the interpretation of future CMB data is 
not limited by knowledge on Ho- 

For a nearby spiral galaxy at distance D (in Mpc) 
which is inclined by i degrees, with a rotation speed of V c 
kms -1 , the proper motion due to rotation equals fine = 7^75 
iiasyr -1 , where k « 4.74 [kms _1 ]/[AU yr -1 ]. For M33 
[i ~ 56°, D ~ 0.84,V C ~ 97], M31 [i ~ 77°, D ~ 0.77, K ~ 
270] and the LMC [i ~ 35°, D ~ 0.055, V c ~ 50] we find; 
H%3 3 ~ 24, nMS 1 ~ 74 and /x^ c ~ 192 /iasyr" 1 . Thus, 
the rotational motions of these galaxies are easily resolved 
by SIM, and, depending on magnitude, also by GAIA and 
In addition to the absolute size of the rotation speed, 
the relative contribution of random motions of the stellar 
population (<r) are important. The Vc/o values of 27, 9.7 
and 2.5 for M31, M33 and the LMC are likely to be indica- 
tive of the ease with which accurate RP distances can be 
determined. 

The "principal axes" (or Mm) variant of the RP 
method applies to galaxies dominated by circular rotation. 
Individual stars are identified along the major (M) and mi- 
nor (m) axes at similar distances from the centre (with the 
same rotational velocity. Proper motions (/x) on the minor 
axis (Mm) measure the circular velocity divided by D. On 
the major axis, /i^, equals projected by cosi. The cor- 
rection for the systemic motion (indicated by superscript 
"C") can be approximately achieved when two stars are cho- 
sen per principle axis, one on each side of the centre (e.g., 

I l^M I = I MM, approaching f^M ,receeding \ Radial Veloc- 

ities (V r ) on the major axis give the cir cular velocity pro- 
jecte d by sini. Combining these yields (|Peterson fc Shaol 
Il997h : 



COS I = 



Mm 



1 c 



D* 



(Mm) 2 



(17) 
(18) 
(19) 



However, the principal axes method requires stars close to 
the principal axes, making it difficult to find enough targets, 
while the effects of major perturbations such as warps, spiral 
arms etc. are difficult to handle. 



9 We performed extensive simulations to convert accuracies in 
position or parallax to proper motion errors, with as a variable 
the inverse of the mission duration (xt). We find that e M /e p0 s ~ 
1.901 x t + 0.474 af, while e M /e w Fa 2.751x4 + 0.981a;?. Thus, for 
mission durations of 5 years, a position error at V=18 of 5.38 
tias for SIM correspond to proper motion errors of 2.1 /xasyr - 1 . 
SIM's best grid accuracy is 3.47 /j,as at V=ll, so that 1.35 



collecting area, quantum efficiency and integration time, that it 
can achieve parallax Eros of 26, 45 and 78 tias at V=15, 16 and 
17, respectively. This leads to proper motion errors of 15, 26 and 
46 tiasyr - 1 , at the same magnitudes. 
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This method can be generalized to a star arbitrarily 
positioned in the galaxy. We use rectangular (x and y) and 
polar (R and 9) coordinate systems centred on the galaxy 
centre and co-planer with the galaxy disk. The x and y' axes 
along the major and minor axes, respectively, with y' is the 
foreshortened y coordinate. The following elementary rela- 
tions between the coordinates and the various projections of 
the orbital velocity V c apply: 



14 

Vy 
V r 



-snV c sin 9' 
snV c cos O 1 



-V m 



tan ( 



VySini — snVc cos 9' sin i 
V y i V y cos i sqV c cos 9' cos i 

K D k D k D 

Vc _ — soVcSinfl' 

~kD ~ ~kD ' 



(20) 
(21) 
(22) 
(23) 

(24) 



with V x and V y the projections of the V c vector on the x and 
y' axes. The angle 9' is the angle between V c and V y , while 
sn = -1 (+1) for (counter-) clockwise rotation. 

The three unknowns D, i and V c can be recovered from 
the three observables (V r , fix, and fi y i) by solving eqns. (f22)l 
- ([24]). Equations ((25)) - (|27)> below codify the "individual 
star" (IS) method: 

V r I -y'//V (25) 



Dis 



xfi x + y'Hy 



COS lis = 



y /V 
-X/J.X 



Vo,I S = V r 



/V 



Xf_ly 



xfi x + y'yiyi 



(26) 
(27) 



where we assume circular orbits. 

The achievable distance error depends on distance, in- 
clination, position angle and (almost) linearly on observing 
errors (OP2000). Evaluating eqn. (|25|l at various positions 
(x, y') in the target galaxy, and if we assume that the accu- 
racies for the radial velocity observations (ev r ) an d proper 
motions are of order the internal velocity dispersion of the 
tracer population, we obtain a distance error (per star) of 
13% for M31. For individual stars in M33 and the LMC, 
we get errors per star of 28% and 90%. Note that for M 31 
and M33 the random motions of a young population (10 
kms" 1 ) corresponds roughly to the expected proper motion 
accuracy of SIM of 2.1 ^iasyr _1 but is much smaller than 
GAIA's and OBSS's accuracies. For the LMC, the internal 
motions roughly equal the astrometric errors of GAIA and 
OBS$F\ Also note that the distance error of 28% for M 33 is 



10 For M31 and a 
kpc) ~ 2.7 /xasyr 

- For M33 and a 
kpc) ~ 2.5 /xasyr - 1 ■ 

- For the LMC, the 



e HGAIA I 



and 



-10 kms~ J 
-10 kms" 1 

e MS/M ^ £ CG4JA> 

'errors" are dominated by the velocity dis 



and 



<t/(k770 
<t/(k840 



persion (a ~ 20 kms 1 ) of the stars: 
77 /iasyr" 1 > e p 



<t/(k55 kpc) 



qu ite close to the value obt ained from HbO-masers method 
by iBrunthaler et al. ] (|2005h . 

Thus, for RP programs targeted at either M 33 or M 31, 
the main obstacle in achieving a small distance error is the 
smallness of the rotational signal. SIM, with its small errors, 
is thus much better suited for such a project. On the other 
hand, SIMs small errors would be "wasted" on the LMC 
stars because there the error budget is dominated by the 
dispersion of the stellar population. 



4.1 Realistic Rotational Parallaxes 

If not modeled carefully, non-circular motions due to spiral- 
arm streaming motions, perturbations from nearby galaxies, 
a bar, warps, etc. may result in a biased distance determi- 
nation (OP2000). We have already seen that the obtainable 
errors for the idealized case of stellar-velocity-dispersion ef- 
fects only yields errors of 13%, 28% and 90% per star for 
M31, M33 and the LMC, respectively. Since the errors due 
to random motions are stochastic, these errors would go 
down with the inverse of the square-root of the number of 
"stars" in the sample. In addition to the uncertainty due to 
velocity dispersion, spiral-arm streaming motions can reach 
amplitudes of 10% of the rotation velocity, leading to sys- 
tematic distance errors of similar magnitude. For the LMC 
the effects of the bar and tidal interaction are even more im- 
portant than for M 31 and M 33. In order to achieve errors of 
several percent, it will be necessary to correct for any sizable 
deviations from circular motion. OP2000 indicate that such 
can indeed be achieved with the next-generation astrometric 
satellites. 



4-1.1 Rotational Parallax: General Case 

For the general case, we separate the three-dimensional mo- 
tion of a star in four parts: 1) the systemic motion (V sy s) of 
the galaxy, 2) a random velocity V a , 3) the circularly sym- 
metric part of the orbital velocity V c , and 4) the peculiar 
component V v . The idea behind this separation is to group 
physical effects together that cause large-scale correlations 
between stellar motions that are amenable to modeling. The 
systemic and circular terms are easy to model, while the 
random term can only be taken into account in a statisti- 
cal sense. The ease with which the peculiar velocity can be 
modeled depends on its cause. Thus, the total space motion 
of a star in an external galaxy is given by: 



V total = V ays + V orbit + V a , and where 

V orbit = V c + Vp . 



(28) 
(29) 



Each term projects onto the orthogonal directions (x,y,z). 
V c corresponds to the rotation speed V c , with components 
V c ,x = V c sin9 and V c , y = V c cos#, and tan(9 = y/x. If the 
systemic motion makes an angle i s with respect to the sky, 
the three observable velocities are: 



kD^x = V x = V m 



+ V a . 



+ 



© 0000 RAS, MNRAS 000, 000-000 



Accurate Extra- Galactic Distances and Dark Energy 13 



K.D flyl 



V c , x + V p>x (30) 
= V sys>ry / cosi s + (V a ,z + Vp, x ) sini + 

(V C ,y + Vp,y + Vcr,y) COS 1 (31) 



■ sin i s — {V a ,z + Vp, z ) cos i + 



V., 



(V c ,y + Vp,y + V a ,y) Sill « , 



(32) 



where V avary > is the component of the systemic motion 
in the plane spanned by the radial and minor axis direc- 



the velocity-gradient effects [Gould I (|2000h ; and § 12.1.1) are 
taken into account self-consistently. These relations involve 
three observables (fi x , fi y i and V r ) and 3x3+5=14 unknowns 
(3 Vsys terms, 3 V a projections, 3 V p terms, V c , radius R, the 
azimuthal angle 8, distance and inclination). Furthermore, 
all observables are expressed with respect to a coordinate 
system specified by four additional unknowns: galaxy centre 
(xo,yo), position angle of the major axis (A) and height (z) 
above the midplane of the target galaxM 11 !. Two more ob- 
servables are available in the position (x, y') of the star. We 
end up with 5 observables and 18 unknownV\ 

Obviously, this set of equations cannot be solved. How- 
ever, substantial simplifications can be made when "averag- 
ing" over several stars. First, the velocity dispersion terms 
need not be determined for individual star but only for the 
group as a whole, so that only 15 unknowns remain (3 un- 
knowns shared by all objects). Second, all stars share the 
same systemic motion, so that the V 3ys terms follow from 
"averaging" suitably located objects (3 unknowns shared by 
all objects). Third, stars at similar distance from the galaxy 
centre share the same rotation speed, inclination, origin and 
orientation of the coordinate system (5 unknowns shared by 
all objects) , and fourth, the distance is approximately con- 
stant 13 ! (1 unknown shared by all objects). Equations (|30[) 
(I32p can still not be solved because they involve five ob- 
servables, six star-dependent unknowns (V p , R, 8, and z) and 
twelve shared unknowns (xq, yo, cj>, Va-, V sys , V c , D, and i). 

To make progress, additional assumptions are needed 
to reduce the number of star-based unknowns. Several ap- 
proaches are possible. One can assume that z and V p lie in 
the plane of the circular component as was done by OP2000. 
This assumption is justified if the perturbations arise from 
within the galaxy itself (e.g., from the bar or spiral arms). 
Thus, Vp, z = = z, so that observations of Na stars yield 



11 Although 2 does not appear explicitly in equations (1 30 1 ) - 
(1326 . it is present through the relation between the observable 
position y' and y: y' = y cos i + z sin i 

12 We could add three more unknowns to handle any radial 
gradients in inclination, position angle and rotation speed. 

13 The nearby RP target galaxies are so close that they extend a 
substantial depth along the line of sight: ~7%, ~3.3%, and ~2%, 
for the LMC, M31, and M33, respectively. This problem is solved 
by replacing D by the sum of the distance to the system centre 
(Dsys) and a position dependent distance: D = D syB + d(x, y, z). 
Note that this procedure does not add new unknowns, just more 
complexity. 



4A?a + 12 unknown values and 5Na observables. The assump- 
tions can be generalized by employing simple parameteri- 
zations for V p and the average mid-plane location < z > 
as a function of azimuth. For example, the galaxy may ex- 
hibit somB_comigation of the plane [as in the Milky Way; 
e.g., iLevine et al.l ([2006) , and references therein] in which 
case z may be expressed as a low-order Fourier series with 
Nz,f coefficients, which behave as shared variables. In the 
general case with Nsv shared variables and four star-based 
variables, N* stars yield 4iV* + Nsv unknowns and 5iV* ob- 
servables, leading to the requirement that: 



N* > N sv ■ 



(33) 



Thus, observations of Nsv stars would suffice to determine 
four star-based parameters for each individual star, as well 
as the Nsv shared unknowns. Observing more objects will 
decrease the errors of the to-be determined parameters. Note 
that the equations of condition are mildly non-linear due 
to the trigonometric terms. However, good initial estimates 
are available for all shared unknowns, so that an unbiased 
solution should be attainable via iteration. This solution 
sc heme is more robus t than the simple method employed 
by lOlling fc Peterson I ||2000D. and will be discussed in more 
detail elsewhere ( Oiling fc Peterson llin preparation] ). 

Employing the solution method outlined above, four of 
the six phase-space coordinates can be determined for each 
observed star, while two are modelled. Such knowledge for 
a large number of stars is likely to be sufficient for the de- 
termination of non-circular motions, as well as the galaxy 
distance. As noted in the sections on extra-galactic water 
masers (§§ 12.1.21 and §§ I2.1.3p . the water-maser method is 
limited by the fact that it samples only two lines of sight 
through the galaxy. On the other hand, the rotational par- 
allax method is designed to map the galactic velocity field 
over a large range in azimuth with the explicit goal to be 
able to map out the non-circular motions. As a result, the 
RP method will yield more robust distances than the water- 
maser method, or, as a matter of fact, any other proposed 
method. 



4-1.2 Rotational Parallax: Observational Requirements 

The distance errors can be reduced by observing more stars, 
if systematic errors in the data allow for such averaging. 
However, because non-circular motions can be correlated on 
large scales, a substantial number of stars needs to be used 
to be able to identify and correct for those systematic non- 
circular motions. Furthermore, these stars must be spread 
out over an area that exceeds the region affected by the 
non-circular motions. OP2000 envisaged using stars spread 
around an annulus in the target galaxy and estimated that a 
minimum of 200 stars are required to achieve a 1% distance 
error for M31. Thus, with SIM proper motions for about 
200 stars per system, and if the non-circular motions can 
be handled properly, the distances to M31, M33 and the 
LMC can be determined to 0.92, 2.0 and 6.4 percent. The 
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required SIM observing time for a rotational-parallax pro- 
gram depends on the magnitude of the stellar targets. For 
GAIA, with a fixed integration time per star, the final accu- 
racy is determined by the total number of stars per target 
galaxy, and the achieved astrometric accuracy as a function 
of magnitude. 

If background galaxies and/or QSOs and stars internal 
to the Local Group galaxies can be measured simultaneously 
in the magnitude range 20 - 24 or so, a pointed, wide-field 
instrument such as OBSS or TPF-C might be well-suited to 
obtain the required astrometry. However, for a 1% distance 
estimate for M31 or M33, the imager needs to be stable 
to approximately ^i ga i axy /lQQ/S.t ~ 0.76A£ ^tas, or about 
3.8 /xas for a five year mission (At = 5). This amounts to 
about l/ll,000 i/l of a pixel. The HST experience [1/100" 1 
of a pixel] i ndicates that such accur acies would be hard to 
accomplish l|Anderson fc King 2003|). Thus, an interferome- 
ter such as SIM is the natural instrument to use when very 
high accuracies are required. On the other hand, the large 
angular extend of the LG galaxies ensures that there are 
tens of thousands of usable background galaxies and many 
more galaxy-member stars, so that wide-field imagers might 
be able to the job. 

By comparing on-galaxy star counts with off-galaxy 
counts on fields with the same Galactic latitude, we esti- 
mate that there are a sufficient number of potential targets 
in these three local group ga laxies. Based on the UCAC2 
catalog (IZacharias et aZj|2004 ). we estimate that the LMC 
contains at least 23,000 stars brighter than V=16 within 
two degrees from the cen tre. Similarly, the 2MASS cata- 
log jSkrutskie et al. Il200d ) yields 2,009 (± 265) and 984 (± 
197) stars with K s < 15 for M 31 and M 33, respectively. The 
LMC is a natural target for GAIA-based RP studies because, 
at the limiting magnitude of its radial velocity instrument 
(V ~ 17), GAIA's proper motion accuracy is smaller than 
the LMC's rotational signature of 190 /^asyr -1 . 

iBrunthaler et all (|2006h estimate the number of water 
masers in M 31 to equal about 17 at an SKA sensitivity 
of 1 mjy. Applying their formalism, we expect a similar 
number water masers in M 33. As we will show below, such 
small numbers of H2O masers spread out over the disk of 
a galaxy are probably not very useful for the determina- 
tion of unbiased distances at the percent level. Given the 
likely absence of tidal interactions in M 33, we can proba- 
bly neglect V p , z and z. Because these masers will be spread 
out over the disk, the gradients parameters (for rotation 
speed, inclination and position angle) may also be impor- 
tant, so that the M33 case may require 15 shared variables, 
so that a minimum of 15 water masers are required [cf. equa- 
tion (|33p ]. In such situations one might resort to using the 
resul ts from the analysis of the ne utral hydrogen (H I) emis- 
sion l|Corbelli fc Schneider I [l997h . However, when employ- 
ing those rotation-, inclination- and position-angle curves, 
the solution is no longer self-consistent. Furtherm ore, one 
may n eed to apply the velocity-gradient corrections [Gould I 
(2000); and § I2.1.1| to interpret the HI data. In contrast, the 



large number of bright stars (~1000) makes M33 a natural 
target for SIM-based RP studies. 

It has long been suspected that the interactions of 
M32 and NGC205 cause spiral structure and warpin g 
in M31 \e.q iBvrd I (|l978l . Il983h : ISofue fc Katol <|l98lft ; 
ISato &: Sawa I (| 19861 )]. Recent Spitzer data indicates that 
the so-called "ring of fire" (where most star formation oc- 
curs) is not centred on M 31, and that the observed "split" in 
this ring is l ikely caused by a rec ent passage of M 32 through 
M 31's disk (|Gordon et aZ. IbOOfJ ). Given that most of M 31's 
brightest stars (and hence water masers) are likely to reside 
in the ring of fire, and that the large-scale dynamics of these 
stars/masers may be significantly affected by tidal interac- 
tion, M 31 may be a less suitable target for an RP program. 
Furthermore, the (expected) small number of water masers 
is probably not adequate for a reliable distance measure at 
the percent level. On the other hand, its large angular size 
makes M 31 a natural candidate for stable, space-based wide- 
field imagers such as OBSS. 

A successful dynamical model of th e LMC would require 
at leas t all the ingredients uncovered bv lvan der Marel et al.l 
(2002), which include: 1) a time-dependent inclination, 2) 
an elliptical disk, 3) low rotation speed, 4) large internal ve- 
locity dispersion, 5) and a significant bar. In addition, the 
time-dependency of the position angle of the line of nodes 
needs to be considered. This model is thought to be suffi- 
cient to characterize the tidal effects of the Milky Way on 
the LMC. As such, the LMC may be a less favorable target 
for a rotational parallax program. Neglecting these pertur- 
bations, the minimum number of required stars (for a 1% 
distance) to overcome just the velocity-dispersion effects ex- 
ceeds: N^lmc > (90%/l%) 2 = 8, 100. However, given the 
much larger number of available bright stars, an accurate 
rotational parallax for the LMC might still be achievable 
with GAIA. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we discussed a new, extremely accurate 
rung of the extra-galactic distance ladder, namely the 
rotational-parallax method. Our preliminary analysis in- 
dicates that the best results might be obtained for the 
following galaxy/instrument(s) combinations: M33/ SIM, 
M 31 1 SIM I OBSS and LMC/GAIA/OBSS. In particular, we 
expect that the M 33-SIM combination can obtain a bias- 
free distance with errors no larger than a percent or so. 
For M31/ SIM the attainable accuracies are likely to be 
dominated by the degree to which the tidal distortions can 
be modeled. The results for galaxy/ OBSS-combinations are 
limited by the level of systematics of the focal plane array. 
For the LMC/ GAIA the main uncertainty is the degree by 
which localized errors can be beat-down by \ // ~N effects. 

Assuming y/~N statistics for the Rotational Parallax 
method, 2% distance errors may be achieved with SIM-like 
observations of just 43 stars in M31 and 200 stars in M33. 
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The LMC requires at least 2,025 stars to reach the desired 
distance accuracy, mainly because random motions are sig- 
nificant, so that it is a preferred target for survey missions 
such as GAIA. 

An accurate determination of the "first" rung of the 
extra-galactic distance ladder does not yield the Hubble con- 
stant. Distance indicators that extend well beyond the Local 
Group must be used for that. Thus, the RP method would 
enable a bias-free calibration of the secondary calibrators. 
On the other hand , the Water Maser Cosmology project 
l|Braatz et al\ \2006) might be better than the RP method 
in the short term because it can be already applied to a 
number of galaxies, albeit with less accuracy per system. 

We expect that GAIA's delivery of five-dimensional 
phase-space coordinates for a magnitude limited sample of 
individual stars of external galaxies (LMC, SMC et cetera) 
will have many significant galactic-dynamics applications 
such as bar-dynamics, disc dynamics and tidal interactions. 
On the other hand, M31 and M33 naturally lead to SIM- 
based programs because the internal motions (relative to 
the stellar velocity dispersions) are four to ten times smaller 
than for the LMC. 

In this paper we focussed our cosmological analyses on 
the case that PLANCK data will be the norm for CMB 
data, and other data (Baryon Acoustic Oscillation, Galaxy 
Clusters, Supernovae and Weak Lensing) have accuracies 
at a level coined Sta ge I by the Dark Energy Task Force 
l|Albrecht et alj|2006h . In such a situation, the effects of de- 
creasing the error on the Hubble constant on our knowl- 
edge of the equation of state of Dark Energy (and Qtot) 
declines when Ho is known more accurately. For example, 
with PLANCK-like CMB data, decreasing the error on Ho 
by a factor of two from the current uncertainty will almost 
half the error on w. However, when decreasing ejj from 2% 
to 1%, e w (and £n tot ) decreases by only 20% (see Fig. [TJ. 

Scaling our analysis to the situation that the other data 
set are much more accurate than at present, we find that 
the effects of smaller en values become gradually less im- 
portant. In accordance with the DETF we find that when 
the quality and quantity of the other data sets reaches the 
Stage IV level, decreasing tH has only a minor effect. 

An alternate utility of determining Ho with great accu- 
racy from nearby galaxies is that it is fully independent of 
traditional cosmological analyses which determine Ho, while 
employing only a weak prior of its value. Especially when 
both methods yield Ho values with accuracies at the percent 
level or so, it becomes possible to test fundamental assump- 
tions underlying CBM-type analyses such as the validity of 
general re lativity as a descript ion for the expansion of the 
Universe l|Albrecht et alj|2006h . 

Other fundamental test of cosmology will become avail- 
able via SIM, GAIA and OBSS astrometry in the form of 
determinations of very accurate ages of stars in (eclipsing) 
binary systems (jLebreton 1120001 ). For such systems, all rel- 
evant fundamental stellar parameters (mass, radius, lumi- 
nosity and metallicity) can be determined with great ac- 



curacy, leaving no wiggle room for stellar models. The ac- 
curacy of stellar dating is fundamentally limited by lumi- 
nosity (distance) uncertainty because stellar evolution pre- 
dicts luminosity (L) evolution of about AL/L « (0.24 — 
0.106 log tms) ± 0.05 per Gyr, where tms is th e main- 
sequen ce (MS) lifetime of the star in Gyr [e.g., lOlling I 
(2003)]. A G6V star (0.89 M ) has a MS lifetime roughly 
equal to the age of the Universe, while its luminosity evo- 
lution amounts to about +12% per Gyr. If one-percent dis- 
tances (2% in L) are available for such stars in the Galactic 
halo, then a measurement of their apparent magnitude cor- 
responds to a determination of their age to approximately 
1000 x (2%/12%) = 166 Myr. Such binary star observations 
will also yield the stars' interior Helium abundance (Y), so 
that Y(t) can be determined and extrapolated to t = 0. 
Many thousands of such halo stars are accessible by GAIA, 
SIM and OBSS, so that a very accurate lower limit of the 
age of the Universe can be obtained right here in the So- 
lar neighbourhood. If one-percent distances are available for 
nearby Local Group systems, such accurate age determina- 
tions can potentially be made in those galaxies, so that a 
wide range in star-formation histories can be sampled. 

I thank Dean Peterson, David Spergel, Ed Shaya and Alan 
Peel for useful discussions. I am also grateful for several sug- 
gestions by the referee to improve the paper, especially with 
regard to Hu's paper and the need to distinguish between 
an accurate extra-galactic distance and an accurate Hubble 
constant. 



© 0000 RAS, MNRAS 000, 000-000 



16 Rob P. Ollmg 



REFERENCES 

Albrecht, A., et al. 2006, |astro-ph/0609591| 

Alves, D. R. 2004, NewAR, 48, 659 

Alves, D. R., 2000, ApJ, 539, 732 

Anderson, J., & King, I. R. 2003, PASR 115, 113 

Argon, A. L., Greenhill, L. J., Moran, J. M., Reid, M. J., Menten, K. M., & Inoue, M. 2004, ApJ, 615, 702 
Armstrong, J. T., et al. 1992, AJ, 104, 241 

Avelino, R P., Martins, C. J. A. P., Nunes, N. J., & Olive, K. A. 2006, PhRvD, 74, 083508 
Bellazzini, M., Ferraro, F. R., Sollima, A., Pancino, E., & Origlia, L. 2004, A&A, 424, 199 
Bellazzini, M., Ferraro, F. R., & Pancino, E. 2001, ApJ, 556, 635 
Benedict, G. F., et al. 2007, AJ, 133, 1810 
Benedict, G. F., et al. 2002, AJ, 124, 1695 

Blake, C. A., Abdalla, F. B., Bridle, S. L., & Rawlings, S. 2004, New Astronomy Review, 48, 1063 
Bom, F. R., Sparks, W. B., & Macchetto, F. D. 1999, A&AS, 138, 253 
Bonanos, A. Z., et al. 2006, Ap&SS, 304, 207 

Braatz, J., Greenhill, L., Condon, J., Henkel, Ch., Lo, F., Reid, M., 2006, the "Water Maser Cosmology Project" at: 

http://www.cfa.harvard.edu/wmcp/index.html and 

http://www.cfa.harvard.edu/wmcp/talks/braatz.legacy.mtg.pdf 
Brunthaler, A., Henkel, C, de Blok, W. J. G., Reid, M. J., Greenhill, L. J., & Falcke, H. 2006, A&A, 457, 109 
Brunthaler, A., Reid, M. J., Falcke, H., Greenhill, L. J., & Henkel, C. 2005, Science, 307, 1440 
Byrd, G.G., 1983, ApJ, 264, 464 
Byrd, G.G., 1978, ApJ, 226, 70 

Carroll, S. M. 2001, Living Reviews in Relativity, 4, 1 

Cooray, A., & Seto, N. 2005, ApJ, 623, L113 

Corbelli E. & Schneider S.E., 1997, ApJ, 479, 244 

Da Costa, G. S., & Armandroff, T. E. 1990, AJ, 100, 162 

Davis, J., et al. 2005, MNRAS, 356, 1362 

Di Benedetto, G. P. 2002, AJ, 124, 1213 

Draine, B. T., & Bond, N. A. 2004, ApJ, 617, 987 

Durrell, P. R., Harris, W. E., & Pritchet, C. J. 2001, AJ, 121, 2557 

Edberg, S.J., Shao, M., Beichman, Ch.A. (editors), 

SIM PlanetQuest, A Mission for Astrophysics and Planet-Finding 

http://planetquest.jpl.nasa.gov/documents/WhitePaper05verl8_final.pdf 

or from: http:/ /planetquest.jpl. nasa.gov/SIM/sim_science_goals.cfm 
Efstathiou G., Lawrence C, Tauber J. and the Planck Science Team, 2005, "The Planck Blue Book," 

www.rssd.esa.int/SA/PLANCK/docs/Bluebook-ESA-SCI(2005)l.pdf 
Eracleous, M., Livio, M., Halpern, J. P., & Storchi-Bergmann, T. 1995, ApJ, 438, 610 

ESA, 1997, The Hipparcos and Tycho Catalogues, ESA Publication Division, ESTEC, Noordwijk, the Netherlands 
Feast, M., Catchpole, R.M., 1997, MNRAS, 286, LI 

Flambaum, V. V. 2006, AIP Conf. Proc. 869: "Atomic Physics 20," p. 29 

Fitzpatrick, E. L., Ribas, I., Guinan, E. F., Maloney, F. P., & Claret, A. 2003, ApJ, 587, 685 

Fomalont, E. B. 2005, ASP Conf. Ser. 338: Astrometry in the Age of the Next Generation of Large Telescopes, 338, 335 
Fomalont, E., & Reid, M. 2004, New Astronomy Review, 48, 1473 
Freedman, W. L., et al. 2001, ApJ, 553, 47 

Gezari, S., Halpern, J. P., & Eracleous, M. 2007, ApJS, 169, 167 

Gordon, K. D., et al. 2006, ApJ, 638, L87 

Gould, A., & Uza, O. 1998, ApJ, 494, 118 

Gould, A. 2000, ApJ, 528, 156 

Gould, A., 1995, ApJ, 452, 189 

Gould, A., 1994, ApJ, 426, 542 

Groenewegen, M. A. T. 2000, A&A, 363, 901 

Groenewegen, M.A.T., & Oudmaijer, R.D., 2000, A&A, 356, 849 

Groenewegen, M. A. T., Romaniello, M., Primas, F., Mottini, M., 2004, A&A, 420, 655 
Guinan, E. F., et al. 1998, ApJ, 509, L21 

Herrnstein, J. R., Moran, J. M., Greenhill, L. J., & Trotter, A. S. 2005, ApJ, 629, 719 



© 0000 RAS, MNRAS 000, 000-000 



Accurate Extra- Galactic Distances and Dark Energy 17 



Herrnstein, J. R., et al. 1999, Nature, 400, 539 

Herrnstein, J. R., Greenhill, L. J., & Moran, J. M. 1996, ApJ, 468, L17 
Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 
Hu, W. 2005, ASP Conf. Ser. 339: Observing Dark Energy, 339, 215 

Humphreys, E. M. L., Argon, A. L., Greenhill, L. J., Moran, J. M., & Reid, M. J. 2005, ASP Conf. Ser. 340: Future Directions 

in High Resolution Astronomy, 340, 466 
Johnston, K. J., et al. 2006, PASP, 118, 1428 

Kaluzny, J., Stanek, K. Z., Krockenberger, M., Sasselov, D. D., Tonry, J. L., & Mateo, M. 1998, AJ, 115, 1016 
Karachentsev, I. D., et al. 2006, AJ, 131, 1361 

Kervella, P., Bersier, D., Mourard, D., Nardetto, N., & Coude du Foresto, V. 2004, A&A, 423, 327 
Kervella, P., Nardetto, N., Bersier, D., Mourard, D., & Coude du Foresto, V. 2004, A&A, 416, 941 
Sakai, S., Ferrarese, L., Kennicutt, R. C, Jr., & Saha, A. 2004, ApJ, 608, 42 

Kurucz, R. L. 2002, ASSL Vol. 274: New Quests in Stellar Astrophysics: the Link Between Stars and Cosmology, 3 
Lane, B. F., Creech-Eakman, M. J., & Nordgren, T. E. 2002, ApJ, 573, 330 
Lebreton, Y. 2000, ARA&A, 38, 35 

Lee, M. C, Freedman, W. L., & Madore, B. F. 1993, ApJ, 417, 553 
Levine, E. S., Blitz, L., & Heiles, C. 2006, Science, 312, 1773 
Lewis, A., & Bridle, S. 2002, PhRvD, 66, 103511 

Macri, L. M., Stanek, K. Z., Bersier, D., Greenhill, L., & Reid, M. 2006, 2006, ApJ, 652, 1133 
Madore, B.F., & Freedman, W.L., 1998, ApJ, 492, 110 
Madore, B. F., & Freedman, W. L. 1995, AJ, 109, 1645 

Makarov, D., Makarova, L., Rizzi, L., Tully, R. B., Dolphin, A. E., Sakai, S., & Shaya, E. J. 2006, AJ, 132, 2729 

Mitchell, R. C, Baron, E., Branch, D., Hauschildt, P. H., Nugent, P. E., Lundqvist, P., Blinnikov, S., & Pun, C. S. J. 2002, 

ApJ, 574, 293 
Mouhcine, M. 2006, ApJ, 652, 277 

Mouhcine, M., Rich, R. M., Ferguson, H. C, Brown, T. M., & Smith, T. E. 2005a, ApJ, 633, 828 
Mouhcine, M., Ferguson, H. C, Rich, R. M., Brown, T. M., & Smith, T. E. 2005b, ApJ, 633, 810 

Newman, J. A., Ferrarese, L., Stetson, P. B., Maoz, E., Zepf, S. E., Davis, M., Freedman, W. L., & Madore, B. F. 2001, ApJ, 
553, 562 

Ngeow, C, & Kanbur, S.M., 2006, ApJ, 642, L29 

Nordgren, T. E., Lane, B. F., Hindsley, R. B., & Kervella, P. 2002, AJ, 123, 3380 
Nugent, P., et al. 2006, ApJ, 645, 841 
Oiling R.P., Peterson D.M., in preparation 

Oiling, 2003, "Connecting the Physics of Stars, Galaxies and the Universe," 

http://www.astro.umd.edu/~olling/FAME/Astrometry_Stars_Galaxies_Universe.pdf 
Oiling R.P., Peterson P.M., 2000, |astro-ph/0005484| 

Paczyhski, B., 1997, in "The Extragalactic Distance Scale," ed. M. Livio, M. Donahue, & N. Panagia (Cambridge: Cambridge 
Univ. Press), p. 273 

Paczyhski, B., & Sasselov, D. 1997, "Variables Stars and the Astrophysical Returns of the Microlensing Surveys;" Ed. by R. 

Ferlet, J-P. Maillard and B.e Raban, Cedex, France: Editions Frontieres, 1997., p. 309 
Panagia, N., Gilmozzi, R., Macchetto, F., Adorf, H.-M., & Kirshner, R. P. 1991, ApJ, 380, L23 

Panagia, N. 1999, IAU Symp. 190: Edited by Y.-H. Chu, N. Suntzeff, J. Hesser, & D. Bohlender. "New Views of the Magellanic 
Clouds," p. 549 

Peebles, P. J., & Ratra, B. 2003, Reviews of Modern Physics, 75, 559 
Perlmutter, S. 2005, Physica Scripta Volume T, 117, 17 
Perryman, M.A.C. et al. 2002, ApSS, 280, 1 

Peterson, D. & Shao, M. 1997, ESA SP-402: Hipparcos - Venice '97, 402, 749 
Iben, I., Jr., & Renzini, A. 1983, ARA&A, 21, 271 

Rizzi, L., Bresolin, F., Kudritzki, R.-P, Gieren, W., & Pietrzyhski, G. 2006, ApJ, 638, 766 

Romaniello, M., Patat, F., Panagia, N., Sparks, W. B., Gilmozzi, R., & Spyromilio, J. 2005, ApJ, 629, 250 

Reid, I. N. 1999, ARA&A, 37, 191 

Reid, M. J. 1993, ARA&A, 31, 345 

Ribas, I., Jordi, C, Vilardell, F., Fitzptrick, E.L., Hilditch, R.W., Guinan, E. F. 2005, ApJ, 635, L37 
Ribas, I., Jordi, C, Vilardell, F., Gimenez, A., & Guinan, E. F. 2004, New Astronomy Review, 48, 755 
Riess, A. G., et al. 2004, ApJ, 607, 665 

Ribas, I., Fitzpatrick, E. L., Maloney, F. P., Guinan, E. F., & Udalski, A. 2002, ApJ, 574, 771 
© 0000 RAS, MNRAS 000, 000-000 



18 Rob P. Oiling 



Salaris, M., & Cassisi, S. 1997, MNRAS, 289, 406 

Sandage, A., Tammann, G. A., Saha, A., Reindl, B., Macchetto, F. D., & Panagia, N. 2006, ApJ, 653, 843 

Sasselov, D.D., et al., 1997, A&A, 324, 471 

Sato, N.R., Sawa, T., 1986, PASP, 38, 63 

Schilizzi, R. T. 2004, Proceedings of the SPIE, 5489, 62 

Seljak, U., McDonald, P., & Makarov, A. 2003, MNRAS, 342, L79 



The SIM book, 1999 http://sim.jpl.nasa.gov/library/book.html 



Skrutskie, M. F., et al. 2006, AJ, 131, 1163 
Soderblom, D. R., Nelan, E., Benedict, G. F., McArthur, B., Ramirez, I., Spiesman, W., & Jones, B. F. 2005, AJ, 129, 1616 
Sofue, Y. & Kato, T., 1981, PASP, 33, 449 
Sparks, W. B. 1996, ApJ, 470, 195 
Sparks, W. B. 1994, ApJ, 433, 19 
Spergel, D.N., et al, 2003, ApJS, 148, 175 
Spergel, D.N., et a/., 2006, astro-ph, 063449 

Stanek, K. Z., Kaluzny, J., Krockenberger, M., Sasselov, D. D., Tonry, J. L., & Mateo, M. 1998, AJ, 115, 1894 
Storchi-Bergmann, T., et al. 2003, ApJ, 598, 956 
Strateva, I. V., et al. 2003, AJ, 126, 1720 

van der Marel, R. P., Alves, D. R., Hardy, E., & SuntzefT, N. B. 2002, AJ, 124, 2639 
White, M. 2006, New Astronomy Review, 50, 938 
Wyithe, J. S. B., & Wilson, R. E. 2002, ApJ, 571, 293 
Wilson, R.E., 2004, NAR, 48, 695 

Zacharias, N., Urban,S. E., Zacharias, M. I., Wycoff, G. L., Hall, D. M., Monet, D. G., & Rafferty, T. J. 2004, AJ, 127, 3043 



© 0000 RAS, MNRAS 000, 000-000 



Accurate Extra- Galactic Distances and Dark Energy 19 



APPENDIX A: WATER MASERS & 
ELLIPTICAL ORBITS 

For the general case of elliptical orbits, Newtonian gravity 
and a dominant central object, the orbital speed at distance 
R is given by: 



V (R) = 



G M 2 - a 



V - km / s (R) = 2ttk 



M 2 - a 



a :AU 



(Al) 



(A2) 



where G is the gravitational constant, M the mass of 
the central object, a the semi-major axis of the orbit, 
and a = R/a . The subscripts "c," "e," and "o" indicate 
cases that are valid for circular, elliptical and general or- 
bits, respectively. We also specify the units of the parame- 
ter in the subscript. For eqn. (|A2f) we express the velocity 
in kms -1 , with the mass of the central object always in 
units of Mq, k (~4.74) converts velocities in AU yr _1 to 
kms -1 and a -Au equals a in astronomical units. For cir- 
cular orbits (a = 1), eqn. (|A1[I reduces to the well-known 
V (R) = <J G M/a , while for orbits with eccentricity e at 
we have V ; P eri = yG M/a a x (l + e)/(l — e) at peri-centre 

and V ;apo = \/G M/a x (1 — e)/(l + e) at apo-centre. 

For our test geometry (see §§ 12. 1 .3|> . the high- velocity 
masers are located at the t angent of the ellipse at distance 
Rhv = \Jb 2 + (ea e ) 2 = y/ (1 — e 2 ) al + (ea e ) 2 = a e from 
the centre, where b is the semi-minor axis of the ellipse. 
Thus, for the high-velocity masers, a equals unity, so that 
their observed radial velocities map out a rotation curve that 
is identical to the case of circular orbits. 

To derive the expected kinematics, we use a coordinate 
system with x and y' axes that are aligned with the apparent 
major and minor axes (neglecting the position-angle warp), 
and with y the corresponding coordinate in the plane of the 
circum-nuclear disk (i.e., y = ycosi). Because the maser 
spots move only a small amount in azimuth during a multi- 
year experiment, we can approximate the orbital motion as 
being locally sinusoidal with constant angular velocity, and 
with period P ~ 2 -nR/V . Then, for Heliocentric distance D, 
the positions, radial velocities (V r ), radial acceleration (V r ) 
and proper motion (fi) are: 



V r ,HV;o 
Smi H V;o 

| Vr,3;o | 

sinis ;o 

I Vr,S;o I 
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R sin 



R cos ( 
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Ro sin 
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Rs,o 
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(A3) 
(A4) 
(A5) 

(A6) 

(A7) 

(A8) 



where we use the additional subscripts "S" and "HV" to 
indicate whether the expression is valid for the systemic or 
high-velocity masers, respectively. The units of the proper 
motion are masyr -1 if the velocities are in kms" 1 and the 
distance in kpc. 

First, we will examine the constraints that derive from 
the high-velocity maser spots. For our toy model, the veloc- 
ities in eqns. (|A2[I and (|A5f) are equal so that the observable 
radial velocity is: 



V r . 



— 1-kk sin iHV o 



Mo 1000_D o 



;Mpc 



(^HV;o;m.as 



(A9) 



where aHV:o;maa is the observed major-axis position of the 
high-velocity spot. Bringing all the observables to the left- 
hand side (LHS) yields: 

aHV;o;masVr HV]0 .2- „ n / A 1 n\ 

— — = Sin lHV;oM D -M V c ■ ( A 10) 

40007T 2 K 2 

Turning our attention to the systemic masers, we use the 
observed minor-axis position of the systemic maser spots to 
determine an expression for the semi-major axis. First, we 
have: y' s . o;AU = Rs-,o-,au cos i s - = a . A u(^ - e D )cosi S;o = 

y's;o;mas X (1000 D a ; Mpo) , SO that: 



CLS;o;AU 



Q>S ;o;mas 
y S :o;mas 



yS;o;mas X ( 1000 -Do; A/pc) 

(1 — e ) cos %S\o 

1 



(1 — e D ) cosis-,o 



(All) 
(A12) 



where the second relation has all observables on the LHS. 
Inserting eqn. (|A11|I into ()A2|) and the result into (|A8|) . we 
find: 



10^ , 

^^_2 yS;o\mas X ^S;o;mas/yr 



M cosis ;o (l + e D ) 

TP 



(A13) 



Further information is contained in the observed accelera- 
tion of the systemic maser spots: inserting eqn. (|A2|I evalu- 
ated at peri-centre into eqn. (|A7|) and using (|A11[) we obtain: 



Vr ; g;o(1000i/s ;O ;, 
47TK 2 



sin i s -o cos 2 i s - M a (1 + e p ) 
LP 



(A14) 



If one assumes circular orbits, a constant inclination and 
that as ~ a,HV, then there are two sets of three equations 
[[ (|AT0)) . ((ATI)) and (|A"T3)) ] and [ (|ATo)) . (fAT2)l and (|AT4)) ]] 
that can be used to solve for the three unknowns: distance, 
inclination and the mass of the central object. If a simple 
warped model [i = io + (3(R — Ro)] is adopted, one can de- 
termine the slope (/3) from the high-velocity maser spots be- 
cause they have a warpe d (projected) geometry and a small 
non-Keplerian behavior (jHerrnstein et al~lll996l ). 

Equations fATO)) , (|AT2|) . (fAT3)> and (fAT^ can be rear- 
ranged into the following relations between the circular and 
eccentric cases: 



sini H V;cy/ M c D C ;M P c = sin iHV-.e \f M e D e - Mp c (A15) 
cosis ;c = (1 — e)cosis ;e (A16) 
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M c cosis-,c M e cosis ;e (l + e) ,. . 



D z D 3 

c;Alpc e;Mpc 

sinis;c cos 2 is ; c M c sinis ;e cos 2 is ;e M e (1 + e) 
Efl - IP 

c;Mpc e;Mpc 



(A18) 



where constraint equations (|A15[l - (|A18j) also equal the LH 
sides of equations (|A10|) . (|A12|) . (|A13[) and (|A14I) . respec- 
tively. Again we can include a fifth relation from the warp- 
ing behavior, but for illustrative purposes we will assume a 
constant inclination, so that we have four equations [ (|A15I) 
- (|A18|) ] with four unknowns [D, i, M and e]. Expressing the 
first three unknowns as a function of the fourth, we find: 

^£ = e + 2e 2 + 0(e 3 ) (A19) 

cos i c 1 — e 



t = V(r^" 1+2e+ 3 e2+0<eS) (A20) 



f " V(^~ 1+4e+ ? cI+0(c3) ' (A21) 

where we note that these relations hold for systemic masers 
at peri-centre. One should keep in mind that eqns. (|A19[) 
- (|A21|) were derived for the special geometry described in 
§§ 12.1.31 These simplifications were used to illustrate the ef- 
fects of unmodeled eccentric orbits for nuclear water masers, 
not to provide accurate relations for any specific case. 

Above we have worked out the case that the systemic 
masing occurs at peri-centre. For the alternate case (masing 
at apo-centre), we find that we can replace e with — e in 
equations (|A11|I - (|A2 If) . including the approximations of 

(EM - 
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